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ABSTRACT 

We compute burst spectra and time structures arising from synchrotron and inverse Compton 
scattering by non-thermal electrons accelerated in shocks which form during the interaction 
between a thin ultra-relativistic fireball and a stationary external medium. We investigate 
the effect of varying the most important model parameters on the resulting burst spectra, 
and we present a set of correlations among the spectral and temporal features of the bursts. 
The spectral hardness, various spectral-temporal correlations and the spectral evolution of the 
simulated bursts are compared to those of observed bursts for a representative set of model 
parameters. Multi-pulse structures are simulated using a variable magnetic field and anisotropic 
emission, and the most important spectral and temporal properties of these pulses are compared 
with observations. 



Subject headings: gamma-rays: bursts - methods: numerical - radiation mechanisms: 
non-thermal 



Introduction 



The observed isotropy (Meegan et al. 1992, Briggs et al. 1996), inhomogeneity (as shown by (V/V max )) 
and the deviation from a —1.5 slope power-law of the logiV — logP distribution for the fainter bursts 
( Meegan et al. 1992] , Horack & Emslie 1994) provide string support for the hypothesis that Gamma- Ray 
Bursts (GRBs) are of cosmological origin. Other observations, such as spectral hardness-brightness 



et al. 1994 



correlations ( Mitrofanov et al. 1992 , Paciesas et al. 1992 ), spectral hardness-duration anti-correlation 
(Kouvcliotou et al. 1993), and the possible time dilation and duration-brightness anti-correlation ( |Norris 



Norris et al. 1995 ; see however Mitrofanov et al. 1996 ), while more equivocal, are also generally 
The large energy that the cosmological source must release suggests 
In this paper we consider bursts that arise when 



compatible with this hypothesis, 
that relativistic effects are likely involved in GRBs. 
an ultra-relativistic cold shell (fireball) is decelerated by interaction with the interstellar medium, or a 
pre-ejected slow wind. As a result of the deceleration an ultra-relativistic blast wave ("forward shock" - 
FS) propagates into the external medium (EM), transferring a substantial part of the fireball kinetic energy 
to the shocked EM, while another shock ("reverse shock" - RS) propagates back into the fireball. This 
is the generic model usually referred to as the "external shock model" ( Meszaros fc Rees 1993| ). If the 
shell is in the linear broadening regime before it is substantially decelerated by the EM (as described by 
Meszaros, Laguna & Rees 1993), a situation that is expected under a wide range of conditions, then the RS 
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is quasi-newtonian and therefore less efficient than the blast wave in converting the shell's kinetic energy 



into heat. In another very likely scenario (the "internal shock model" - Rees Sz Meszaros 1994, Paczyhski 



fc Xu 1994 ) the energy conversion takes place when several ultra-relativistic shells collide with each other, 



before the deceleration caused by the EM becomes important. Here we focus on the first model. 

In order to simulate the propagation of the two shocks and to model the fireball-EM interaction 
at large Lorentz factors (Tq > 100), we have developed a one-dimensional hybrid (finite differencing 



exact Riemann solver) hydrodynamic code ( |Wen, Panaitescu fc Laguna 1996| ). As the conversion of 



kinetic to internal energy takes place, the heat stored in the post-shock gas can be released as radiation, 



generating a burst. In a previous paper (Panaitescu et al. 1997) we have simulated burst light-curves 



from fireballs with moderate Lorentz factors (ro < 200), using a simplified prescription for the energy 
release. The results were single-hump bolometric light-curves with a large temporal asymmetry (light-curve 
decay lasting substantially longer than its rise), practically insensitive to variations of the EM density. 
In order to carry out an appropriate comparison of this model with the rich observational database that 
has been accumulated by BATSE and other experiments, we need to compute the spectra of such bursts, 
and to study the burst spectral evolution and its correlation with the other observational properties and 
parameters of the model, by taking into consideration specific energy release mechanisms. The spectral 
hardness of the observed bursts and its time evolution are well studied, and to reproduce these should 
represent a major goal of any GRB model. In this work, we calculate the effect of the source evolution on 
the burst spectrum, and explore the spectral-temporal correlations predicted by the model. We also explore 
the physical requirements necessary in order for this model to produce multiple-humped light curves, 
and discuss the possible limitations. Spherical symmetry is assumed for simplicity throughout the paper, 
which also describes well the case of jets with an opening angle 9 > T^ 1 . The importance of a non-planar 
symmetry can be assessed from the shape of the light-curves and pulses presented below. 



2. Model Parameters, Assumptions, Approximations, Scaling Relations 

The most important parameters that describe the dynamics of the fireball-EM interaction and the 
energy release mechanism are listed in Table 1, together with most relevant equations in which they appear. 
The evolution of an impulsive fireball has two phases: a free expansion phase, when the amount of swept up 
EM is small and the deceleration caused by it can be neglected; and a decelerated expansion phase, when 
the fireball kinetic energy is used to heat the swept-up EM. The fireball dynamics during the first stage was 



calculated analytically and simulated numerically by Meszaros et al. (1993) . The evolution during this stage 



is determined by three parameters: (1) the energy Eq = 10 51 £q,51 ergs deposited in the ejected fireball; (2) 
the entrained baryonic mass M, parameterized through the dimensionless entropy Tq = Eq/Mc 2 3> 1; and 
(3) the initial size of the fireball ro (which may be of the order of the neutron star radius , ro ^ 10 6 cm). At 
the beginning of the free expansion phase the fireball is accelerated as the radiation energy Eq contained in 
it is adiabatically transformed into bulk motion energy, and becomes stretched out into a thin shell. The 
absence of a strong burst precursor with a quasi-thermal spectrum suggests that most of the initial internal 
fireball energy is in the form of baryonic kinetic energy when the shell becomes optically thin and photons 
escape from it. Therefore, the maximum Lorentz factor attained by the fireball is ^ Tq, corresponding to a 
kinetic energy ^ Eq. 

After the fireball Lorentz factor attains its maximum value the fireball coasts at constant Tq and later 
reaches optical thinness. The deceleration caused by the interaction with the EM must be taken into 
account when the energy stored into the shocked EM is a substantial fraction of the initial kinetic energy 
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TABLE 1 

Summary of the most important parameters and physical quantities that characterize the 
dynamics of the interaction fireball-em and the burst energy release 



Symbol 


Definition 


Equation 


fro 


initial fireball kinetic energy 




r 


initial fireball Lorentz factor 


E = T Mc 2 


tdec 


hydrodynamic time-scale 


© 


Xb 


magnetic field parameter 


Ub = ^sUint 


B 


magnetic field intensity 


i) 


hi 


electron-to-proton energy ratio 


i) 


7m 


minimum electron Lorentz factor 


i) 


£ SY/IC 


co-moving SY/1C photon energy 


®>© 


E SY/W 


detector SY/1C photon energy 




yRS/FS 


Kompaneets parameter 






synchrotron cooling time-scale 





Eq ■ The shocked EM internal energy is much larger than its rest mass energy, since its random (or thermal) 
Lorentz factor is ~ Fo 3> 1. Throughout most of this paper we assume that the EM is homogeneous, 
characterized by a single parameter: its number density n = 1 uq cm~ 3 . The deceleration time-scale in the 
stationary frame (with respect to the Earth) of the center of explosion (the laboratory frame) is 

tdec = r dec /c ~ {Eo/Tlnmpc 5 ) 1 / 3 ~ 8.3 x 10 5 E 1 ^ n 1/3 F" 2 / 3 s , (1) 

where Lq = 10 2 Lq^- Due to the relativistic motion of the source, the stationary observer receives radiation 



emitted in dt in a much shorter time dT — dt/[2T (t)] ( Rees 1966 ), where T(t) < Tq is the Lorentz factor 



of the shocked emitting medium. The burst duration is then approximately 

T b w 10 t dec /2T 2 = 420 El% n~ 1/3 L^ /3 s , (2) 

where a factor of 10 was included in order to account for the progressive decrease of the flow Lorentz factor 
of the radiating medium. Equation (^) and the observed GRB durations imply that 100 Tq few x 1000. 
It also shows that the burst peak flux F p satisfies F p (~ EqD^ 2 ^ 1 ) cx E^ 3 T^ 3 n 1 / 3 D~ 2 , where D is the 
distance to source, if most of the available energy Eq is radiated. 



The dynamics and energetics of the deceleration phase were calculated by Rees & Meszaros (1992) and 



by 3ari & Piran (1995). For computational efficiency, the numerical simulations presented here were started 
from 0.5 tdec, when only ~ 12% of the EM mass within 1 rd ec had been swept up and the deceleration prior 
to this time can be safely neglected. At t = 0.5 td ec the only physical parameter that depends on r is the 
internal pressure P of the fireball, and in fact this pressure is irrelevant as long as the shell is cold (P <C pc 2 , 
p is rest mass density). Therefore, the hydrodynamics of the shell-EM collision is characterized by the set 
of three parameters (Eo,To,n). 
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In the co-moving frame, the shocked EM has typical densities ~ 10 3 particles/cm 3 , and can radiate 
away its internal energy through synchrotron radiation (SY) and inverse Compton (IC) scattering of the 



SY photons, in the presence of a modest magnetic field. Such mechanisms were considered by Meszaros, 



Rees fc Papathanassiou (1994) , who studied the spectral properties of bursts arising from external shocks, 
and by Sari, Narayan & Piran (1996), who derived constraints on the radiation mechanisms parameters 
from the variability observed in most bursts and from efficiency considerations. The galactic magnetic field, 
even when compressed behind the FS, would be too weak to lead to efficient radiation. However, a frozen-in 
magnetic field present in the fireball (and thus in the fluid behind the RS) would allow the post-FS material 
to cool by IC scattering of the SY photons coming from the post-RS medium. The swept-up EM could 
radiate even more efficient if a random turbulent magnetic field builds up in it. In our calculation we use 
for simplicity this latter scenario; moreover, a frozen-in magnetic field will usually have only a fraction of 
the strength of a turbulent magnetic field at equipartition (when the magnetic field energy density is equal 
to the internal energy density of the gas) . 

We use the following assumptions and approximations in order to simulate the emission of SY and IC 
photons from the gas behind the two shocks: 

[1] the magnetic field B is parameterized relative to the internal energy density Ui n t : Ub = As Ut n t , 
where Ub — B 2 / 8ir is the magnetic field energy density. For strong shocks, equation (8) derived by 



Blandford fe McKee (1976) yields Ui n t = 3 x 10 3 tiqT 2 fs erg/cm 3 , where Tps is the Lorentz factor of 



1/2 1 /2 

the FS shock, so that B = 0.27 A^ n Q < T FS G. Since the post-shock fluids are very close to hydrostatic 
equilibrium, Ui n t, and therefore B, have almost the same values behind both shocks. 

[2] shock acceleration leads to a power-law distribution of electrons 

dA4(7e) = Cle V d le , 7m < 7e < 7M , (3) 

where j e is the random electron Lorentz factor and 7V e is the number density of electrons. Such a 
distribution is initialized in every grid cell after it is swept up by one of the two shocks, and its subsequent 
evolution is determined solely by the SY and IC losses. Therefore we do not take into account adiabatic 
losses or further energy exchange between protons and electrons. The former simplification is justified by 
the fact that the electron cooling time-scale is much lower than the dynamic time-scale (as shown below). 
We have taken 7m/7to = 10 because the cooling time-scales for larger 7m would be too short (Ss 10~ 5 ^ ec , 
typically) and it would require a large computational effort to follow accurately the evolution of these very 
energetic electrons. The effect of a larger ratio 7m /7m on the burst spectrum can be easily estimated in the 
figures presented in the next section. Moreover, if To is not low (ps 100) or if the magnetic field is not weak 
or the shock acceleration inefficient (i.e. low j m ), the most energetic electrons radiate at energies above the 
upper limit of the BATSE window (10 keV - few MeV). The power-law index p was chosen to be 3. 

[3] The minimum electron random Lorentz factor7 m is determined by a parameter k which is the ratio 
of the energy in electrons and that in the monoenergetic protons (with Lorentz factor 7 p ), after shock 
acceleration: 

/ dAf e {~f e )m e (7 e - 1) = Kn p m p (j p - 1) , (4) 

An, 

where n p is the density of protons. Using the equality of the sum of the electronic and protonic partial 
pressures and the total pressure (determined by the hydrodynamics of the fireball-EM interaction), j m and 
the electron distribution are completely determined: 

lm {n) = 3 [(p - 2)/(p - 1)] [(1 - X l ~ p )/ (1 - X 2 -?)} f K (m p /m e ) (P/pc 2 ) , (5) 
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where X = Jm/Jth an d /« = «/(«+ 1). This result is valid for 7 P ^> 1 (equivalent with P ^> pc 2 , 
which is true for the fluid behind the FS); a similar result can be obtained in the limit 7 P - 1 < 1 (i.e. 
P <C pc 2 , which is correct for the fluid behind the RS). For the ultra-relativistic FS, equations (8)-(10) 



from |Blandford fc McKee (1976)| lead to P/pc 2 = 0.24 T FS and "f m ,FS = 660/ K r f5 . [Meszaros et al. 
(1993) have shown that the evolution of the fireball thickness A during the free expansion phase and for 
r > ro Tq = 10 10 7-0,6 Tq 2 cm is A = r / Tq. This leads to a fireball density at r ~ rd ec that is much larger 
than that of the EM and produces a mildly relativistic RS. Numerically, we found that the Lorentz factor 
of the RS in the frame of the yet un-shocked fluid is practically independent of Tq: Trs — 1.1 • In this case 
it can be shown that P/pc 2 ~ 4 x 10 -2 which leads to jm.RS — 100 f K . 

[4] In the co-moving frame, the SY radiation emitted by any electron is approximated as monochromatic, 
with a frequency equal to the peak frequency v c (averaged over the pitch angle) of the SY spectrum emitted 
by an electron with Lorentz factor 7 e (t) (that evolves in time, as the electron loses energy): 

e SY ^ 4.0 x 10" 9 7 2 (t)B eV . (6) 

The electron cooling and continuous electron injection will produce spectra that are flatter than the 
spectrum of the SY radiation emitted by a single electron below and above the peak frequency v c 
(vF v (x 1s 4 / 3 and vF v cx v z / 2 exp [— f/z/ c ], respectively), so that the effect of integrating over time and over 
electron distribution hinders the features of a single electron spectrum. Thus, this approximation is in fact 
better than it seems at first sight. 

[5] The spectrum of the SY photons up-scattered in the Thomson regime is also approximated as 
monochromatic, at the average energy of the IC spectrum for 7 e (^) : 

eic = 4/3 7 2 (i)esy ■ (7) 

The Klein-Nishina (K-N) effect on the scattering of SY photons with energies comparable or larger than 
m e c 2 /7e(t) is taken into account. The SY energy density Usy, necessary for calculating the IC losses, 
is computed as an integral over the volume of the shocked media of the SY local output. There is a 
strong relativistic beaming of the SY photons due to the radial motions of the origin of a given photon 
and the place where the scattering takes place: as seen from the co-moving frame of the up-scattering 
region, the SY source is moving away, unless the two regions (of SY emission and of IC scattering) are 
moving in the same radial direction. We assumed that the Usy spectrum is monochromatic, at the peak 
frequency of the SY spectrum generated by the most numerous (and least energetic) electrons that are 
in the same volume element where the IC scattering takes place. This approximation is justified to some 
extent by the aforementioned strong relativistic beaming and the geometrical dilution of the SY output, 
which should make the contribution to the Usy of the SY emission from the vicinity of IC scattering place 
to be dominant. Due to this assumption the IC spectra shown in the next section are calculated using 
only the following combinations: (i) SY-RS photons scattered on electrons accelerated by the RS and (ii) 
SY-FS photons scattered on electrons accelerated by the FS. The mixed combinations (iii) SY-RS photons 
scattered by FS electrons and (iv) SY-FS photons scattered by RS electrons are not taken into account. We 
will assess the effect on the computed spectra of neglecting the last two combinations. 

The approximation of taking the SY and IC spectra of a single electron and that of the SY photon 
field to be up-scattered as monochromatic is done for computational efficiency. The results presented in the 
next section were obtained from numerical runs that last few hours on a Sparc Sun 20 Station for the lowest 
r considered there (To — 100) and up to few days for the highest T we used (To — 800). Most of this 
computational effort is used to calculate the burst spectrum by integrating over the volume of the shocked 
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fluid (which reduces to integrating over the radial coordinate and the angle relative to the line of sight 
toward the center of symmetry) and over the electron distribution in each infinitesimal volume element, 
and repeating this triple integral after a time short enough to accurately treat the evolution of the most 
energetic electrons (which have the shortest cooling time-scale). Adding another integral within the triple 
integral, in order to include the real SY or IC spectrum from a single electron, would lead to excessively 
long runs. 

The above analytic considerations and approximations allow us to calculate the energy E p at the peak 
of the power per logarithmic energy interval (vF u ) for the SY and IC spectra from both shocks, as seen 
from the detector frame. Numerically, we found that about 50% of the total energy released by a burst 
is emitted from t = lt c i ec until t = l.btdec- During this time Tps decreases from ~ 0.6 Tq to ~ 0.4 To, 
therefore, to a good approximation, Tps — To/2, so that 

~/ m ,p S ~ 1-3 x 10 5 /« (r /400) (8) 

and 

B ~ 54 X^ 2 n /2 (r /400) G . (9) 

Taking into account that the relativistic motion of the radiating fluid boosts the co-moving energy by a 
factor between T (if the fluid moves at an angle T _1 from the line of sight toward the fireball's center) and 
2T (if the fluid moves on this line of sight), where T ~ 0.7Tps is the flow Lorentz factor of the shocked 
fluid, we obtain for the SY-RS radiation: 

E SY,RS „ QA f 2 A V2 n l/2 (rV400)2 ^ (1Q) 

The co-moving energy of the SY-RS photons Ssy,rs = 2 x 10~ 3 f 2 A]/ 2 nj 2 (r /400) eV is well below the 
limit for K-N scattering m e (?/"fm,RS = 5/^T 1 keV, therefore 

Ef > RS ~ 6 ft X 1 ' 2 nl /2 (r /400) 2 keV . (11) 

The SY photons emitted by post-FS electrons with Lorentz factor 7 TOj fs arrive at detector at 

E^ FS ~ 800 f K X 1 ' 2 n X J 2 (r /400) 4 keV , (12) 

and, in the co-moving frame, are too energetic to be up-scattered in the Thomson regime: 
£sy,fs = 4/2 X 1 ^ 2 nl' 2 (r /400) 3 keV » m e c 2 / 7m = 4/- 1 (ro/400)" 1 eV, as long as 

logr ,2 + |log/ K + ilogA B + ilogno > -0.15. (13) 

Let's assume that inequality (^3|) is satisfied (we argue below that it must be so for an efficient burst). 
The energy m e c 2 /j e is a good measure of the photon energy above which the K-N reduction is very 
effective, in the sense that it drastically reduces the intensity of the IC component. The up-scattered 
radiation will be emitted when the FS electrons have cooled enough so that £sY,Fs(je) < m e c 2 /j e , 
implying 7 e < jkn = 5 x lO 4 ^ -1 / 3 , where B is lower than estimated above (eq.||), as the shocked 
material has lost some internal energy. To a good approximation, this fraction can be taken 1/2, so that 
lKN = 1.5 x 10 4 A B 1/6 n~ 1/6 (ro/400)^ 1 / 3 . The SY radiation emitted by FS electrons cold enough to 
scatter their own SY photons in a mild K-N regime have a detector frame energy less than 

C S ^^« n^ 6 (Fo/400) 4 / 3 keV , (14) 
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which gives the peak energy of the up-scattered spectrum from the FS: 

E^ C > FS « 0.6 \ B 1/6 n- 1/6 (r /400) 2 / 3 TeV . (15) 

The optical depth for Thompson scattering of the shocked fireball is trs — E gtk / 47r m p c 2 r 2 ec T — 
10~ 6 El^ x nj 3 (ro/400) 1 / 3 < 1. The effect of IC scattering on the SY-RS spectrum and on electron cooling 
can be assessed through the Kompaneets parameter Y RS = Jm.RS T RS ■ 

y b.s „ io- 2 /^ 1 / 5 3 1 n^ 3 (r /400) 1 /3 « i . (16) 

Calculating a similar Kompaneets parameter for the FS is more difficult because earlier accelerated 
electrons can be so cold that they scatter their own SY photons in the Thomson regime while more recently 
accelerated electrons are very energetic and scatter their SY photons in the extreme K-N regime. A 
simple way of obtaining upper limits for this parameter would be to assume that all electrons have the 
same random Lorentz factor and that the up-scattering takes place at the limit between the Thomson and 
the K-N regimes. (At given energy So of the incident photon, the Kompaneets parameter for electrons 
colder than 70 = m e c 2 /eq increases as 7 2 while for electrons with random Lorentz factors above 70 the 
same parameter increases as ln[27 e /7o].) The SY photons emitted by electrons with 7 TOi fs (eq. are 
up-scattered in this mild K-N regime by electrons that have 7 e = 140 f~ 2 A^ 1 ^ 2 n Q (r /400)~ 3 ; for such 
scatterings the Kompaneets parameter is: 

yFS „ 1Q - 5/ -4 A -l £l/3 n -l/a (ro/400) - 2 0/3 <<; 1 

Before reaching 70 calculated above, electrons are cold enough to scatter the SY photons they produce (see 
equation for jkn above). Assuming again that all electrons are monoenergetic and have 7 e = Jkn, the 
Kompaneets parameter is: 

Y^^lO-^-^E^n^iTo/AOO)-^- (18) 
For electrons colder than jkn the Y parameter should increase as 7 2 while for more energetic electrons the 



same parameter should decrease as 7,T 4 . Thus equation (18) gives an upper limit on the expected intensity 
of the IC-FS component relative to that of the SY-FS emission. In deriving equations (|l^) and (|l^) we 
approximated the mass of the swept up EM by a fraction 1/Tq of the fireball mass. We can conclude from 
equations ([l6]) and (|T^) that the IC emission is not expected to alter substantially the intensity of the SY 
radiation from the two shocks or the synchrotron cooling time-scale of electrons. 

The energy release mechanisms considered in this model involve only two important parameters (k,\b) 
which, based on equation (112), must satisfy the double inequality 



0.1 < lo g r , 2 + ~ log/ K + i log A B + i log n < 0.7 , (19) 

to ensure that the burst fluence in the BATSE window corresponds to a significant fraction of the total 
energy radiated by the source. In the laboratory frame, the synchrotron cooling time of the least energetic 
FS electrons is t SY = (3 m e c / 4 a Th U B j m ) T ~ 140 f- 1 A^ 1 n^ 1 ( r /400 )~ 2 s, or 

t SY ^Ax 10- 4 f- 1 X B ' n- 2/3 E-K 3 (ro/400)- 4 / 3 t dec . (20) 

If the SY cooling time-scale is larger than the hydrodynamic time-scale tdec, the progressive fluid deceleration 
and the adiabatic cooling of the shocked fluid lead to a softening of the spectrum (less energetic electrons + 
lower Doppler blueshift) and reduce the burst intensity (less energy radiated away by electrons). The end 
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result is a weak, soft and possibly un-detectable burst. Therefore, efficiency considerations also require that 
t SY < tdec, which, using equations ([!]) and (pcf), leads to 

3 If 

logr , 2 + -log(/ K A s ) + -log^o.si + 2 logn ° > ~ L9 ■ ( 21 ) 

Note that if this condition is satisfied by the electrons with the minimum random "f m then it is also satisfied 
by the more energetic electrons (with j e > 7 m ). 

The observed burst durations determine the range of To (from cq. || and 10 ms ^ T), ^> fOOOs, it results 
that 100 £ To < 5000), thus equations ( JTgj ) and (|2f]) can be used to constrain the energy release parameters 
(As, ft). Numerical simulations for fireball Lorentz factors Tq > fOOO require a large computational effort, 
so hereafter we will restrict our attention to cases with To < 1000, which give burst durations ^ f s 
(from eq. 0]), i.e. those bursts that are most often considered in the GRB statistics. For such initial 
Lorentz factors, k must be larger than f0 -2 and Xb should not be less than 10~ 4 in order to give a 
spectral peak in the BATSE window. The burst fluence in the BATSE window is determined also by the 
fraction of the available energy Eq which is radiated at a power large enough to give at detector a photon 
flux above a given threshold. It would be wrong to assume that this fraction is strictly proportional to 
f K = ft/(ft + 1), the fractional energy in electrons after shock acceleration because, even if electrons and 
protons are completely "decoupled" after shock acceleration (i.e. no further energy flow from protons to 
the rapidly cooling electrons), the heat stored in protons drives forward the FS, which accelerates new 
electrons. In this indirect way a substantial fraction of proton energy can be transferred to electrons and 
radiated. Numerically we found that in At = 2 tdec a burst with k = 0.1 and Xb = 1 radiates ss 50% of the 
total energy E , which is not much less than the « 80% of E that a burst with k = 1 and A^ = 1 radiates 
during the same time. For this reason it can considered that k does not have an important effect on the 
energy released as long as it is not much less than 10 -1 . 

It is easy to see that if most of the SY-FS radiation is in the BATSE window (i.e. eq. [[l9| is satisfied) 
then either the K-N effect reduces severely the IC emission (eq. [[f3) is fulfilled) or the Kompaneets parameter 
Y FS (eq. @) is less than 1. This means that if a burst observed by BATSE represents the SY radiation 
emitted by the shocked EM, then the IC-FS radiation from the same fluid is less energetic than the SY-FS 
emission and can be safely neglected in calculating the cooling time-scale. On the other hand, if equation 
( |f3| ) is not satisfied (i.e the K-N cut-off does not reduce the efficiency of IC scattering behind the FS) then 
the SY-FS radiation does not arrive in the BATSE window. This suggests that a burst visible to BATSE 
can also be obtained from the IC-FS radiation if the efficiency conditions [1] 10 keV < E^ C ' FS < few MeV, 
[2] Y FS > I and [3] t 1 ^ < td ec are simultaneously satisfied (condition [3] is relevant for the burst efficiency 
only if condition [2] is satisfied). It can be easily shown that condition [1] implies up-scattering in the 
Thomson regime and that it cannot be fulfilled at the same time as condition [2]. In other words, any 
combination of parameters 100 Tq ^ 5000, 10~ 2 ^k<1, As < 1> Eq,5i ~ 1, no ~ f leads to either an 
IC-FS component that contains a substantial fraction of the available energy but is at energies larger than 
those visible to BATSE, or to an IC-FS radiation that arrives mainly in the BATSE window but is much 
less energetic than the SY-FS radiation emitted by the burst, due to a small Kompaneets parameter. Thus 
the observed bursts could be IC-FS radiation only the initial fireball kinetic energy is much larger than 
10 50 ergs/sr, in which case a much more energetic emission should be detected at energies lower than the 
BATSE window. Similar conclusions have been reached previously by Sari et al. (1996). 

In principle two other model parameter constraints can be obtained if it is required that all electrons 
are confined in the shocked fluid and that the duration t acc of the of the electron acceleration process is 
much shorter than the corresponding SY cooling time-scale, ensuring that electrons can reach factors "f e 
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larger than the post-FS j m derived above. The former condition requires the electron gyration radius 
R g = j e m e c / eB to be less than the thickness A (<; rdec/^o) of the shocked fluid shell, while the latter 
condition requires that t acc ~ R g /c <§; t SY . It can be shown that if the inequalities (|l9| ) and are 
satisfied then electrons are indeed confined in the shocked fluid and are accelerated on a time-scale much 
shorter than the SY cooling time-scale, so that these two conditions do not bring any new constraints on 
model parameters. 

The total set of model parameters used in the simulations below is (Eg, Tq, n; k, Xb] D), including the 
luminosity distance to source, to which the parameter 7A//7m could be added in some special cases. In the 
following numerical results D = 10 28 cm; however the cosmological redshift effect was not accounted for 
specifically, because the spectral redshift and temporal dilation can be included in the model independent 
of the hydrodynamic simulation. Therefore, the temporal-spectral correlations discussed below are not of 
cosmological origin; they are intrinsic properties of the bursts. 



3. Results and Comparison with Observational Data 

A hardness-brightness correlation, hardness-duration anti-correlation, and brightness-duration 
anti-correlation are straightforward predictions of these external shock models. From equations @, 
and the fact that the peak flux scales as F p oc E^ 3 T 8 / 3 n 1 / 3 D~ 2 we see that the fireball's initial Lorentz 
factor r (100 < T < few x 10 3 ) is the parameter with the strongest influence on the spectral and temporal 
burst properties. If the other parameters have a relatively narrow range (1 < -Eo,5i ~ 10, no ~ 1) or 

are within the limiting values suggested above (0.1 < K < 1, 10~ 4 < \b < 1), then the correlations or 

3/2 

anti-correlations expected among the burst parameters are due to their To-dependence, and are E p cx F p , 
E p oc T. , and F p oc T, . Evidence for a hardness-brightness correlation has been presented Mitrofanov 



et al. (1992), Paciesas et al. (1992), |Ncmiroff ct al. (1994)| , [Pclaez et al. (1994)| , and |Mallozzi et al. (1995)| , as 



it is implied by the hardness ratios, break energy or E p dependencies on the peak count rate or brightness 
class shown in these articles. A quantitative comparison is not easy as authors seldom use F p and E p 
in their analyses (or at least the same definition of the burst hardness); nevertheless it appears that the 
observed correlation is weaker than predicted above. The hardness-duration anti-correlation is observed by 



Dezalay et al. (1992) and Kouveliotou et al. (1993) (see however Band et al. 1993), while the evidence for a 



brightness-duration anti-correlation is controversial ( Norris et al. 199E , Mitrofanov et al. 1996); if present 



it is probably far weaker than indicated by the above analytic scaling. Of course, a distance dispersion 
of an order of magnitude, as well as a broad luminosity function (variations of Eq, Xb and n parameters 
among bursts) and evolutionary effects would all tend to mask such an F p — T& anti-correlation through the 
parameter Tq. 

Further comparison with observational data can be done using numerical results. Figure 1 shows spectra 
(computed as flux weighted averages of 10 instantaneous spectra, uniformly distributed within Tf,) generated 
with different values of Tq when the other parameters are held constant. The IC component from the RS 
is shown separately while the other components can be distinguished in the spectrum and are identified in 
this figure. Note that most of the burst energy is in the SY component from the FS and that an important 
fraction of this energy arrives at detector in the BATSE window if the parameters Xb and k are close to their 
maximum values (as predicted by eq. [ [l9[ ) . The ratio 7m /j m is relevant for the burst fluence in the BATSE 
window only for Tq = 100. The spectra cover a fairly broad range in energy (13-15 orders of magnitude). 
The burst spectral flux at 550 nm is w 10" 10 (r /400) 8 / 3 ergs cm" 2 s" 1 eV" 1 = 40 (r /400) 8 / 3 mJy, which 
corresponds to a magnitude V ~ 13 — 6.7 log (r /400). 
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We can now estimate the effect of approximation [5] above (mixed RS-FS combinations in the 
IC spectrum are neglected), using the previous equations for the minimum electron Lorentz factor 
behind each shock and equations (|l0|) and (|l2|). The energy (in the laboratory frame) of the SY-FS 
photons that would be up-scattered by post-RS electrons with ^ m .Rs a t the limit between Thomson 
and K-N regimes is E^~ FS = (m e c 2 /j m , RS )T ~ If- 1 (r /400) MeV. Equation @ and Figure 
1 show that there are SY-FS photons less energetic than E^ S ^ FS . Therefore, due to approximation 
[5], a fifth component of the spectrum (SY-FS photons IC scattered in the RS) is neglected. This 
component would have a peak below E§ S N ^ FS ^ n RS ~ 10/ K (r /400) GeV if E§ S N ^ FS < Eg Y < FS 
or at EfS^ FS = E^ FS 7 l <RS ~ 10 f* n)j 2 (r /400) 4 GeV if E§jr FS > E* Y > FS . The energy 
of the SY-RS photons that would be up-scattered by FS electrons with j m .FS m a mild K-N regime 
is Eftf^ FS — (m e c 2 /j m ,Fs) T — 0.8 J" 1 keV. Equation ( [To| ) and Figure 1 show that there are 
SY-RS photons at energies lower than E^f* FS , Approximation [5] does not take into account a 
sixth component of the spectrum (SY-RS photons IC scattered in the FS) that would appear at 
Efc^ FS = E^ YRS 7^ FS = Ef^ FS . It can be shown that the cooling of FS electrons with j m ,FS through 
this kind of scatterings is less efficient than through SY emission. Therefore, the numerical results do not 
take into account the mixed components for IC scattering and under-estimate the burst flux in the lower 
energy part of the IC-FS components shown in Figure 1. Fortunately, the flux in the most important energy 
range (the BATSE window) is very little affected. Otherwise, the intensity of the IC component relative to 
the SY emission from each shock, as shown in Figure 1, is consistent with the previous estimations (eqs. 
@, 0, and @). 

The peak energy E p of the spectra shown in Figure 1 passes through the BATSE window as Tq is 
increased from 100 to 800. As expected, higher Lorentz factors lead to harder spectra (see legend). This can 
be also seen using the hardness ratio HR32, defined as the ratio of counts in the third BATSE channel (100 
keV - 300 keV) to that in the second channel (50 keV - 100 keV): HR 32 (r = 100,7 Af /7 m = 100) = 0.46 
(T b ~ 500s), HR 32 (r = 200) = 0.50 {T b ~ 100s), HR 32 (r = 400) = 0.80 (T b ~ 10s), and 
HR 32 (r = 800) = 0.98 (T b ~ 2 s). Figure 2 shows the SY-RS spectra obtained for a fixed r = 400 and 
combinations of parameters (n; As, n) in which only one parameter is changed relative to the "standard" 
combination (1 cm -3 ; 1, 1), showing the effect produced by each parameter and allowing comparison with the 
spectral peaks given by equation (|l2]). The hardness ratios for the new spectra are HR 32 (0.1; 1, 1) = 0.62, 
HR 32 (1; 0.1, 1) = 0.62, and HR 32 (1; 1, 0.1) = 0.43 . 

The hardness ratio range allowed by the model is less wide than the range of E p generated by the 
values of Tq considered, and it is useful to compare these ratios with those of the observed bursts. The 
HR 32 values above are consistent to those presented by Paciesas et al. (1992) and comparable to those 
found by Nemiroff et al. (1994), Mitrofanov et al. (1996), and Kouveliotou et al. (1993). According to 
this last reference, the average HR 32 is 0.87 for bursts with T b > 2 s. The hardness ratios HR43 of the 
simulated bursts range from 0.2 to 0.4 for 200 < F < 400 (10 s < T b < 100 s) and is ~ 0.6 for T = 800, 
in good agreement with the values calculated by Dezalay et al. (1992). For the bursts shown in Figures 
1 and 2 the ratio HR 3 4 i 2 of the photon fluxes in channels 3+4 (above 100 keV) and in channels 1+2 (25 
keV - 100 keV) is between 0.25 and 0.65, lower than the hardness ratios calculated by Bhat et al. (1994): 
0.3 < HR 3 4 i 2 < 1. It is also important to compare with observations the low and high energy spectral 
indices a and /3 as defined by Band et al. (1993): a is the asymptotic limit of the slope of the photon 
spectrum Ce = dN 1 /dE at arbitrarily low energy (Ce oc E a exp [— E/Eo\), and f3 is the slope of Ce at 
energies higher than the spectral peak E p (Ce oc E 13 ). For the Tq = 400 and Tq = 800 spectra shown in 
Figure 1, the spectral indices a obtained from fits of spectra below E p , using the Band function, are —1.8 
and —1.6, respectively. The high energy spectral indices for the same initial Lorentz factors are —2.9 and 
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—2.8 . These values are consistent with those found by Band et al. (1993): —1.5 < a < and —3 < < — 1. 
The expected analytic value of a is —1.5 (integrated spectrum of SY radiation from cooling electrons) while 
that of j3 is —(1 + p/2) = —2.5 (spectrum of SY radiation from a steady-state distribution of electrons with 
continuous power-law injection). The slightly lower values of the indices obtained numerically are due to the 
continuous deceleration of the FS, leading to a progressive spectral softening through decreasing magnetic 
field, Doppler blueshift factor and random Lorentz factor of injected electrons. The increased steepness 
of the SY-FS spectra shown in Figure 1 below ~ 1 keV is due to IC scattering in Thomson or mild KN 
regimes, as predicted by equation (|l4|). 



A spectral evolution of GRB from hard to soft has been observed by many authors (e.g. Norris et al. 



1986| , |Mitrofanov et al. 1992| , |Band et al. 1992| , [Bhat et al. 1994| , [Ford et al. 1995| ). Figure 3 shows the 



light-curve and temporal evolution of the spectrum resulting from a simulation with constant parameters 
Xb and k. A substantial fraction (60%) of the burst radiation falls within the BATSE channels 1-4. The 
burst light-curve exhibits a sharp rise and a slow decay during which the flux is well approximated as a 
power-law F oc T -12 . The bottom graph shows the burst's hard to soft spectral evolution: the hardness 
ratio HR32, the mean energy E m in the BATSE channels 1-4 (defined as the ratio of the energetic flux and 
photon flux in this band), and the peak energy E p decrease monotonously during the burst (see legend). 
During the light-curve decay (T > 3 s), these spectral parameters can be approximated by power-laws in 
T: HR32 oc T -01 , E m oc T~ 02 and E p oc T~ 12 (similar indices describe the spectral evolution of the other 
bursts shown in Figures 1 and 2). The peak flux and spectrum of this burst show that its peak photon 
flux in the BATSE window is of order 0.1 7/ cm 2 s, corresponding to a weak burst. This is due in part to 
the conservative choice Eq = 10 51 ergs over Air steradians and to the almost maximal luminosity distance 
D = 10 28 cm ~ 10 Gly in this example. Beaming of the fireball in a solid angle < 1 sr would easily boost 
the peak photon flux of this burst above 1 7 / cm 2 s. 

In an efficient burst, the synchrotron cooling time of the FS electrons is much shorter than the 
hydrodynamic time-scale. Consequently, most of the burst radiation is emitted by the leading edge of 
the expanding shell of shocked fluid, from a region which is td ec /t SY ~ 10 3 4- 10 4 times thinner than the 
shell containing all the shocked fluid. At detector time T corresponding to i, the observer is not receiving 



radiation from this very thin sub-shell, but from a very elongated ellipsoid (see Rees 1966) of semi-major 
axis ~ l^dec- Consequently, the detector receives radiation that was emitted at times spread over ~ l^ec, 
which means that the spectrum and light-curve reflect the long time-scale variations of the burst physical 
parameters while all features arising from short time-scale variations are well mixed and less distinguishable. 
The Lorentz factor T(t) of the shocked fluid is monotonously decreasing; therefore, at constant energy 
release parameters n and Xb (relaxation of this assumption is considered in the next section), the spectral 
evolution of the burst shows only the time-evolution of T(t). Thus, the hard to soft spectral evolution 
shown in Figure 3 is purely due to the deceleration of the radiating fluid. 



Burst Substructure 



We further test the ability of the blast wave model to accommodate some of the more frequently 
observed features of spectral evolution in bursts that exhibit individual pulses: 

(i) the spectrum hardens before an intensity spike, and softens while the photon flux is still increasing 
( |Mitrofanov et al. 1992| , |Kouvcliotou et al. 1992] , |Band et al. 1992| , |Bhat et al. 1994j |Ford et al. 1995| ) ; 



(ii) the hardness of successive spikes decreases (Norris et al. 1986 



Band et al. 1992 



(iii) pulses peak earlier in the higher energy bands (Norris et al. 1986 



Ford et al. 1995); 



Kouveliotou et al. 1992 



Norris 
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et al. 1996); 

(iv) pulses exhibit faster rises at higher energies and longer decays at lower energies ( Norris et al. 1996 ) 
and thus peaks are shorter at higher energy (Link. Epstein, fc Pricdhorsky 1993, Fcnimore et al. 1995, 



Mitrofanov ct al. 1996| ), 



although exceptions from these "rules" are not un-common. Since the flow Lorentz factor of the radiating 
shocked fluid is monotonously decreasing, the simple kinematics of this fluid cannot by itself produce 
spectra showing increasing hardness, nor light-curves containing peaks (assuming spherical symmetry), so 
departures from this simplest case need to be considered in order to explain such features. 



4.1. Temporal variability from EM inhomogeneities 

The pulses that are observed in bursts could have some relation to fluctuations in the EM density, 
denser EM blobs leading to a more intense release of energy. In this scenario, the spherical symmetry is 
lost and a 3D hydrodynamic code is required to perform numerical simulations. The following is a purely 
analytical model of the situation. The duration of each pulse is determined by three factors: 

(1) the projection of the shocked inhomogeneity on the line of sight toward the center of explosion, 
determined by: (la) the laboratory frame thickness Sr of the overheated radiating region and (lb) the 
angle 5 9 = R / r subtended by the shocked blob around its position 6 on the spherical cap from which the 
observer receives radiation, where R is the radius of the un-shocked blob, assumed spherical. R should be 
less than the radius r / T ~ r^ec /T ~ 3 x 10 13 n I ' 3 (r'o/400) _5 <'' 3 cm of the visible spherical cap, or else the 
pulse lasts as long as the whole burst, 

(2) the time it takes to sweep up the entire inhomogeneity, 

(3) the laboratory frame duration of the energy release St ~ t SY (eq. Q). 
The contributions of these factors to the pulse duration are: 

(la) AT Sr = 6r/c~2R/r 2 b c, 
(lb) AT se = 26R/c, 

(2) AT R = R (IV 2 + 9 2 ) I c, and 

(3) AT st = st(T- 2 + 2 )/2 7 

where Tf, is the flow Lorentz factor of the shocked blob, which we will approximate by the Lorentz factor T 
of the rest of the shocked EM, although it must be lower because the inhomogeneity is denser. In calculating 
ATsr above we used the fact that in the laboratory frame the shocked blob material is ~ T 2 times denser than 
before the shock, therefore Sr ~ 2 R / T 2 . Since AT$ r and AT R are of the order R/T 2 c and AT$ S ps R/Tc, 
it results that AT$ r , AT R < AT se . Furthermore, if R > 10 10 ^'Ag 1 ^ ^To^OO) -3 cm, then AT St can 
be neglected relative to ATse- Thus, for 1(T 3 K^X^n^iTo/AOO)- 3 AU < R < 1 n 1/3 (r /400)- 5 / 3 AU 
(assumption 1), ATse determines the duration of the pulse. If R is less than the lower limit set above, then 
one has to consider the contribution of the cooling time to the pulse duration. If R is above the upper limit , 
then the pulse duration is comparable to Tb and it would be impossible to have bursts with more than a 
few pulses. 

In order to derive the distribution P(AT) of the durations of pulses in individual bursts, we assume 
that the co-moving photon number spectrum of the radiation emitted by each blob is a power-law 
dNry = Ce~° de (assumption 2) over a range in energies e m in ~ Smax wide enough that the blue-shifted 
corresponding laboratory frame range covers the band in which observations are made, for all blobs that 
are seen by the observer (i.e. for all inhomogeneities that produce at detector a peak photon flux above 
a given threshold Cu m it). Thus we assume that the Doppler shifted edges of the co-moving spectrum: 
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/[r(l - vcos9)] (v = Vr 2 - l/r is the flow velocity) satisfy -E mm < 2? m and 
Em < E max (assumption 3), where E m and £?m are the lower and upper edges of the observational band. If 
so then the peak photon flux C p of each pulse is C p oc [C/(a — l)][r(l — v cos Q)]~ a ~ 2 . The constant C can 
be determined using the fact that the total number of photons emitted per unit time in the co-moving frame 
(= [C/(<7— 1)] s x m ^ nl if a > 1 and e m i n <C s max ) is equal with the number of emitting electrons N e multiplied 
with the number of photons emitted per unit time by each electron, which is independent of the Lorentz 
factor of the electron, and depends only on the magnetic field B. If all blobs are identical not only in size 
but also in density (assumption 4), then N e is the same for all pulses and therefore C/(a — 1) oc Be^. 
The minimum co- moving energy £ m i n of the SY photons is proportional to B and to 7^, where 7 m is the 
minimum Lorentz factor of the electrons accelerated when the FS interacts with the EM inhomogeneity. 
We further assume that the parameters for energy release (k 7 \b,7m /jm,p) are the same for all blobs 
(assumption 5), so that B oc T and C/(a — 1) oc r 3<J ~ 2 . In the end, the peak photon flux at detector for 
any pulse can be written as: 

C P (m= C -^F 1 (E m ,E M ;a)F 2 (p^°) ^ ^ # (1 _ ^g )g+a ■ ( 22 ) 

where F\ and F 2 are generic notations for functions of the indicated variables and nu b is the density of the 
inhomogeneity. 

The condition C P (T,8) > Cumit determines which blobs yield pulses that are detectable, provided 
that the kinematics L(Lo,r) of the shell during the deceleration phase and the spatial distribution 
rib(r) of the blobs are known. For an adiabatic interaction T(r) oc r -3 / 2 ; numerically, we found that 
r(r) = 1/2 To (r/rdec)^ 3 ^ 2 is a good approximation. The number density of the EM inhomogeneities is 
considered to be a power-law: n&(r) oc r~ m (assumption 6), thus the homogeneous distribution is the 
particular case m = 0. Based on these assumptions one can determine for any shell position r the maximum 
angle 9 max (r) relative to the line of sight toward the center of explosion for which C p (r, max [r]) — Cu m it 
and integrate over r and 9 to find the pulse duration distribution. Figure 4 (upper graph) shows this 
distribution for a representative set of parameters (R,To,m,o,C max /Cu m it), where the last parameter is 
a measure of how bright is the pulse from a blob located at (r = 1 rd ec , 9 — 0) relative to the detection 
threshold. The same figure allows one to assess the importance of each parameter: P(AT) is rather 
insensitive to To and depends strongly on R. The lack of correlation with the initial Lorentz factor is 
due to the fact that 9 max (r) is weakly dependent on L while the strong correlation with the blob size is 
clearly implied by ATse oc R. If the observed burst substructure is due to EM inhomogeneities and if the 
assumptions made here are not far from reality, then the latter correlation could be used to infer from 
observations the typical size of these inhomogeneities. As expected, if the pulse detection threshold is 
decreased, longer pulses are seen, as more blobs at larger angles become brighter than Cu m i t ■ 

The statistics of pulses in a set of bursts can be derived by convolving the pulse duration distribution 
for individual bursts with the distribution P(T ) of the initial Lorentz factors of the shells that generated 
the bursts in that set. For this, we assume that P(To) oc Tq for T m i n < T < T max (assumption 7), that all 
shells run into the same EM (assumption 8), and that all bursts distances are the same (assumption 9). 
Such pulse duration distributions are shown in Figure 4 (lower graph) for T min = 200 and T max = 800. 
C m ax (as defined above) for Tq = 200 was chosen 10 times larger than Cumiu this determines C max for 
any other To. It can be seen that P(AT) is not strongly dependent on the parameters v and m. Thus, it 
is possible to estimate the size of the blobs by using durations of pulses in different bursts, as R remains 
the parameter that affects the most the pulse duration distribution. The pulse decomposition performed by 
Norris et al. (1996) shows that for the brightest bursts AT is between 0.1s and few seconds, therefore R 
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must be of order 0.1 AU. 

If radiation is emitted not only from the higher density blobs but also from the rest of the EM, then 
the effect of a more intense emission of radiation from a blob combined with a stronger decrease of the 
flow Lorentz factor induced by the same blobs is likely to lead to a shallow peak. In other words, the 
radiating power of the source is increased but, in the same time, the radiation received by the detector is 
more stretched out in time than the radiation emitted before and after, ironing out the peak. We are forced 
thus to assume that only the blobs emit significant radiation (perhaps due to an enhanced magnetic field). 
In this case, however, a new difficulty arises: as |Sari fc Piran (1997) pointed out, if pulses do not overlap 



significantly then the upper bound on the size of the emitting blobs set by the observed pulse durations 
limits to about 1% the fraction of the area covered by these blobs on the spherical cap visible to the 
observer, leading to a low burst efficiency. (A higher efficiency can be reached if the number of blobs is large 
enough to cover the entire spherical cap visible by the observer, but then the pulses lose their individuality, 
producing a single hump burst.) In order to explain the observed burst fluences, one has then to assume 
that the ejecta is beamed into a fraction 1/100 of the full sky and that almost 99% of the initial energy is 
not released as 7-rays or is lost adiabatically. Thus, in principle, this explanation for a complicated pulse 
structure can work if the ejecta is in a jet, without increasing the total energy above 10 51 ergs, if 99% of 
this energy can go undetected. 



4.2. Temporal variability from energy release fluctuations 

To explore the limits of the ability of external shock models to generate pulses, we consider a second, 
idealized scenario, in which the burst sub-structure is due to fluctuations in the parameters Xb and k which 
characterize the release of the internal energy stored in the shocked gas. Here we consider the case where 
k is constant in time, and we assume a variable magnetic field. A time varying k should have a similar 
effect on the cooling time-scale (t SY oc Xg 1 ft" 1 ), but a stronger one on the spectrum (Ep Y ' RS cx A^ 2 n 2 ). 
If the magnetic fields are such that: (1) at their maximum value the burst radiates mainly in the BATSE 
window and (2) at minimum value, t SY > tdec (the source is in a 7-quiet phase), then equations ( |l2| ) and 
( p0| ) show that Xb must vary by more than 4 orders of magnitude: XB, m in ~ 10~ 4 and XB,max ~ 1; i- e - the 
magnetic field must vary by at least two orders of magnitude. We do not speculate here on the nature of 
the microscopic process that could produce such fluctuations of more than 2 orders of magnitude in the 
magnetic field strength over time-scales that should be shorter than 0.1 tdec, and remark only that plasma 
dynamo mechanisms which build up the field to a fraction of the equipartition value could plausibly result 
in such field variations. In the presence of such variations, multiple peaked bursts are obtained, as shown 
below. 

In Figure 5 we show a burst with two peaks, resulting from a relatively large scale variation of As (see 
the inset of panel a). The spectral evolution is shown with open symbols in panel (b): E m is decreasing 
during the first peak (T Pj i ~ Is), then increases and peaks around T — 2s, approximately Is before the 
photon flux and energetic flux peaks (T p ,2 — 3 s) and monotonously decreases through the remainder of 
the burst. The hardness ratio HR32 shows the same behavior. The monotonous spectral softening of the 
burst during the first peak is due to the deceleration of the shocked fluid and also to the fact that this 
simulation was started from 0.5 tdec- Thus, the radiation emitted by the fluid moving at angles ~ 
(relative to the observer) prior to t = 0.5 tdec is not accounted for, resulting in an artificial softening of 
the spectrum during the first peak that obscures the spectral evolution expected from a variable magnetic 
field. This is not the case with the spectral evolution during second peak, which shows clearly the second 
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As-pulse. The duration and temporal symmetry of each peak can be characterized through the rise and 
fall times T R = J^ p dT f(T), T F = J^ b dT f(T), where f(T) is the photon (or energetic) flux normalized 
to its maximum value (reached at the peak time T p ), and through the time-asymmetry ratio A = Tp/T R . 
T p , the pulse duration AT = Tr + Tp and the ratio A are given for each peak in the legends of panels 
(c) and (d). Note that both pulses are narrower and peak earlier at higher energies, which are features 
known to occur in observed GRBs. The rise and fall times of the pulses decrease with energy, but their 
time-asymmetries show opposite trends: the first pulse appears more symmetric at higher energy while the 
second is more symmetric at lower energies. In log AT — log E, where E is the geometrical mean of the low 
and high edges of the four BATSE channels, the two pulses appear relatively scattered from a straight line; 
nevertheless, if a power-law is fitted, then AT oc T -0 20 . If the pulse full width at half maximum is used, 
then ATpwhm oc T -0 24 . A clearer power-law dependence is found for the single-pulse burst shown in 
Figure 3: AT, ATpwhm oc T~ 015 . Norris et al. (1996) decomposed 41 bright GRBs into pulses and found 
that the average full width half maximum of the pulses varies with energy as E~ 33 if only the separable 
pulses are used, and as E~ 3S for all pulses in the analyzed bursts. Therefore the pulse duration-energy 
anti-correlation of our simulated bursts is somewhat weaker than the observed one. The second peak in 
graph (a) is slightly more time-asymmetric than the first peak (in BATSE channels 1-4: A\ — 5.0 and 
A 2 = 5.4); it also is wider, more shifted to later times at lower energies (graph c vs. graph d) and spectrally 
softer (as shown by HR 32 in graph b). These are exactly the relative features observed by Norris et al. 
(1996) in their pulse decomposition analysis. The blast wave model reproduces the increase in the burst 
hardness before an intensity peak but the simulated spectral hardening is weaker than what is observed. 

If radiation is emitted isotropically in the co-moving frame (as would be the case for a turbulent 
magnetic field), then the observer receives radiation mainly from portions of the fluid moving at angles 
^ r _1 relative to the line of sight. Light emitted by such a spherical cap at time t is spread in detector 
time T over AT(t) = r(t) / 2T 2 (t), where r(t) is the radial coordinate of the cap. Since the flow is 
ultra-relativistic, r(t) ~ ct and thus AT(t) T& (from eq. ||). This means that any instantaneous event 
that occurs in the spherical shell is seen by the observer smeared over a good fraction of the burst. Pulse-like 
emission of radiation and spectral features due to a change in the fluid physical parameters are ironed out 
very efficiently by sampling over the entire opening angle of the region seen by the observer. This naturally 
suggests that, if spherical symmetry in the laboratory frame is maintained, then the angular opening of 
the cap from which the detector receives radiation must be less than T -1 in order to reduce the blending 
of the temporal and spectral features arising from fluctuations in the burst physical parameters. This 
could happen if the radiation, instead of being emitted isotropically in the co-moving frame, is beamed 
along the radial direction of fluid motion. If this radiation is concentrated in two cones of solid angles 
27r(l — /i co ) sr around the radius vector, then the observer receives radiation from a cap of angular opening 

[(i-M C0 )/(i + ^ )] 1/2 r- 1 <r- 1 . 

The effect of such an anisotropic emission can be assessed from Figure 6: as the radiation in the 
co-moving frame is emitted within a narrower solid angle, the light-curve becomes more time-symmetric. 
Due to the monotonous spectral softening (As and k are constant, T decreases), the photon flux decays 
more slowly than the energetic flux and therefore is more time-asymmetric (see the rise and fall times given 
in the legend of each graph). This figure can be compared with similar ones presented by Mitrofanov et 
al. (1996), showing the GRB "average curve of emissivity" in the BATSE channels 2+3. In the isotropic 
case, the radiation emitted by the fluid moving at large angles (< T _1 ) relative to the line of sight is 
Doppler blue-shifted by a factor 2 relative to the radiation emitted by the fluid moving exactly toward 
the observer. This large angle radiation arrives later at the detector and is mixed with the radiation 
emitted at later times, but from regions moving at smaller angles. As the co-moving frame solid angle in 
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which radiation is emitted decreases, the detector receives less radiation from the fluid moving at large 
angles, therefore the radiation emitted at different times is less mixed and the spectrum reflects better the 
instantaneous physical conditions of the radiating fluid. For Figure 6 this means that the spectrum shows 
better the deceleration of the shocked fluid in the anisotropic case than in the isotropic one. This can be seen 
in the evolution of the three spectral parameters used so far during the burst fall (T > T p ): HR32 oc T -01 , 
E m oc T-° 2 and E p oc T" 1 - 2 in the isotropic case, HR 32 oc T~°- 5 , E m oc T~ - 7 and E p oc T~ 2A if in the 
co- moving frame the radiation is emitted within 47r/5 sr around the radial direction while in the most 
anisotropic emission considered here (47r/17 sr around the direction of flow) HR32 oc T~ - 9 , E m oc T~ 13 
and E p oc T~ 2 - 7 . During the burst fall, the Lorentz factor of the leading edge of the expanding gas (from 
where comes most of the radiation received by the detector if t SY <C t^ec) is approximately T oc T~ 2 / 3 , 
which implies that the fastest possible spectral peak evolution is E p (oc T 4 ) oc T~ 27 . Therefore, the most 
anisotropic case considered above yields a spectral evolution that reflects very well the deceleration of the 
shocked fluid. For Figure 5, the anisotropic emission allows the spectrum to show better the effect of a 
varying As (graph b, filled symbols): note the much faster spectral softening during the first pulse and the 
sharp spectral hardening before the second peak. 

The same type of anisotropic emission can be used to generate multi-peak bursts, as shown in Figure 
7. The standard of comparison is an isotropic emission case (top graph of Figure 7), for which the radiation 
coming from the source is blended into a single hump light-curve. The more anisotropic the emission is, 
the shorter and brighter the pulses are and individual peaks can be distinguished better. The progressive 
spectral softening makes these peaks to be less well separated in photon flux than in energy flux, as can be 
seen in the middle graph. Pulses appear more distinct in the case of maximum anisotropy considered here 
(bottom graph). If much internal energy were to accumulate in the shocked fluid between two consecutive 
A^-pulses and if most of it is radiated during a magnetic field pulse, then the observed peaks may be 
blended into a single one (as it happens with the pairs of pulses 1-2, 5-6 and 7-8 in the middle graph). The 
pulse onset times, calculated from the time when t SY < td ec and using the radial coordinate of the shell's 
leading edge, are indicated with numbers. The peak of each pulse occurs slightly later due to the angular 
opening and thickness of the source. Note that later pulses are more time-asymmetric than earlier ones and 
last longer; this is caused by the continuous deceleration of the source. A larger number of pulses can be 
simulated by using an even stronger co-moving frame anisotropy. 



5. Conclusion 

We have discussed some features of numerically simulated GRB spectra and light-curves from external 
shock models, with particular attention to the expected spectral-temporal correlations and the expected 
degree of temporal substructure. The values of the most important model parameters (Fo, n; As, n) were 
chosen such that the burst releases an important fraction of its energy in the BATSE window. No effort 
was made to optimize these parameters so that the simulated bursts mimic the observed ones, other than 
considering (phcnomcnologically) the effects of a variable magnetic field and an anisotropic emission pattern 
in the co-moving frame for some of the models. We then compared the features of the numerical bursts 
with those characteristics of the observed GRBs that are well established, such as the spectral hardnesses, 
low and high energy spectral indices, hard-to-soft spectral evolution, correlation between spectral hardness 
and intensity, pulse duration vs. energy etc. We summarize here the features of the numerically simulated 
model bursts : 

1. The brightness and spectral hardness are correlated. 



- 17 - 



2. They show a spectral hardness - duration anti-correlation : E p oc T b . The observed dependence is 
weaker, which could be due to the variations of the energy release parameters (A#, k) from one burst 
to another, 

3. For single pulse light-curves, the photon flux in the BATSE window rises as T 16 and decays 
approximately as T -10 . The fall is steeper when the co-moving frame emission is anisotropic, 

4. The low energy index of the averaged spectrum is a = —1.7 ± 0.2, not far from the expected value of 
— 1.5 . For p = 3, the high energy index is (3 — —2.8 ± 0.1, not far from —2.5, the theoretical value. 
The former index is determined by the evolution of the accelerated electrons if the spectral peak E p is 
in the BATSE window, while the latter index depends on the choice of the electron power-law index 
P, 

5. The spectra show a general hard to soft evolution outside of intensity pulses. In single hump 
light curves arising from isotropic co-moving frame emission, the spectral evolution at T > T p is 
characterized by E p oc 7 1 - 11 ±0- 1 _ jf m th e co-moving frame radiation is emitted preferentially on the 
radial direction of motion, the spectral evolution is faster, 

6. The peak (or break) energy E p increases with intensity during a pulse, but peaks earlier. The mean 
energy E m in the BATSE window and the hardness ratio HR32 (or similar ones) show a similar trend. 
The increase in the burst hardness before an intensity peak is stronger in the anisotropic emission 
case, 

7. Earlier pulses are harder and have a more time-symmetric profile at higher energies. Later pulses may 
show an opposite trend: more symmetry at lower energies, 

8. Pulses peak earlier and are shorter in higher energy bands than at lower energies. Numerically, we 
found that the pulse duration scales as E~°- 20±0 - 05 , which is a weaker dependence than observed 
(E°- 3 ^ 0A ). We must recognize here that, taking into account some of the approximations made, the 
BATSE channels are relatively narrow for the accuracy of our simulations, so that the calculated 
pulse duration vs. energy dependence can be considered satisfactorily close to what is observed, 

9. The angular opening of the region from which the observer receives radiation limits the number of 
separate pulses to very few. A larger number of pulses results if radiation is not emitted isotropically 
in the co-moving frame. Later pulses are more time-asymmetric than earlier ones and last longer if 
they result from a periodic variation of the source radiating power. 

The above list of burst characteristics is in agreement with, or at least close to, what is observed. It is 
worth noting that the brightness - duration anti-correlation induced by the fireball Lorentz factor Tq will 
be weakened by any dispersion in some of other parameters involved in the model, such as the distance D 
to the burst, the source initial kinetic energy Eq, and the energy release parameters (k, Xb), which could 
explain why this anti-correlation is controversial or, at best, a very weak one. 

Another major observational feature against which to contrast models is the bimodality in duration 



distribution. One reason why this is expected in external shock models of GRB (Sari et al. 1996) is related 
to the limited energy range in which BATSE is sensitive: significant energy arrives at the detector in the 
BATSE window from either the forward shock (FS), in the case of the long bursts, or from the reverse shock 
(RS), in the case of the short bursts, but from neither of these shocks for bursts with durations T ~ 2 s. 
Moreover, such different burst origins can explain the lack of a duration - brightness anti-correlation: the 
RS is less efficient than the FS in converting the fireball kinetic energy into gamma-rays, diminishing the 
brightness of the short burst. In our model, the RS is always mildly relativistic and radiates inefficiently 
(i.e. at energies outside the BATSE window). If the expanding shell thickness increases faster than we 
considered, before its deceleration becomes important, then the density of the colliding shell can be small 
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enough to lead to the formation of a more relativistic RS. The electron Lorentz factor can be further 
boosted if an injection fraction well below unity is assumed. In this case the bursts duration bimodality 



would be reproduced numerically. A different explanation (Meszaros & Rees 1993) for a bimodal duration 
distribution may be that shorter GRB arise from events in a relatively dense external environment (the 
external shocks occur in the progenitor's own pre-ejected wind or in a denser disk galactic disk environment) 
while longer GRB could be due to events in a much lower density environment (e.g. the object has moved 
out of its own pre-ejected wind or it has escaped the galactic disk). 

In summary, the external shock or blast wave model can explain the spectral features and correlations 
of most bursts. It can also explain the time histories of those bursts which have a simple structure (up 
to 4-8 pulses) if the magnetic field is variable and the co-moving emissivity is appreciably anisotropic. It 
is difficult to see how this could be extended to fit also bursts with more than 8-10 pulses. There is no 



difficulty in explaining the latter in outflows with "internal" shocks (e.g. Rees & Meszaros 1994), which are 
expected to have similar spectral properties without limitations on the degree of variability. Nevertheless, 
external shock models show a remarkable degree of qualitative agreement with a large range of medium to 
long time-scale spectral and temporal correlations exhibited by the GRB data. This suggests either that 
external shocks may be responsible for part of the emission of a GRB, or else that a substantial subset of 
bursts (i.e. the less variable ones) may be ascribed to external shock events. 



This research has been supported in part through NASA NAG 5-2362. 
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Fig. 1. — SY and IC spectra for Eq = 10 51 ergs, n = 1cm -3 , D = 10 28 cm, As = 1, k = 1, Jm/j™ = 10, 
p = 3, and different parameters IV The thick dotted curve is for To = 100 and -fAt/jm = 100. Labels 
indicate the origin of each component: SY = synchrotron, IC = inverse Compton scattering, RS = reverse 
shock, FS = forward shock. E Fe — v F v is the power per logarithmic energy (or frequency) interval. Vertical 
dashed lines show the BATSE window. The legend also gives the spectral peak energy E p for each parameter 

r . 
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Fig. 2. — SY spectra from the FS for different sets of parameters (n ; A#,k). _E = 10 51 ergs, r = 400, 
D = 10 28 cm, 7M /7m = 10, andp = 3. E p (l; 1, 1) = 500 keV, £ p (0.1;l,l) = 45kcV, £ p (l; 0.1,1) = 75keV, 
£ p (l; 1,0.1) = 17keV. 
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Fig. 3. — Light-curve (upper graph) and spectral evolution (lower graph) for Eq = 10 51 ergs, n = 1cm -3 , 
To = 400, D = 10 28 cm, Xb — 1, and k = 1. The detector time T is measured from the moment of the 
explosion that generated the fireball. 60% of the burst energy arrives at detector in the BATSE channels 
1-4. The legend of the bottom graph indicates the hardness ratio HR32, the mean energy E m in keV in the 
BATSE window and the spectral peak E p in keV (in this order), at different moments and for the averaged 
spectrum. 
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Fig. 4. — Burst sub-structure from EM inhomogcncitics. Top graph: pulse duration distribution for 
individual bursts. The solid thick curve is for R = 0.1 AU, Tq = 400, a = 2, m = (homogeneous 
distribution of blobs), and a pulse detection threshold Cu m it — C max /W, where C ma x is the photon flux 
that a blob located on the line of sight toward the center of explosion and at r — 1 rd ec yields at detector. 
Other distributions shown are for the same set of parameters except that indicated in the legend. Bottom 
graph: duration distribution for a set of bursts. The combinations of initial Lorentz factor distribution 
and spatial distribution of the EM inhomogeneities are indicated in the legend. Parameters: R = 0.1 AU, 

r m in = 200, T max = 800, (7 — 2, C max {J^min) / C limit = 10. 
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Fig. 5. — Burst sub-structure from energy release fluctuations. Time history, spectral evolution, and pulses 
shapes in BATSE channels 1-4 for E = 10 51 ergs, n = 1cm -3 , T = 400, D = 10 28 cm, k = 1, and time 
variable Xb, shown in the inset of graph (a) (light-curve), (b) Hardness ratio HR32 and mean energy E m in 
the BATSE channels 1-4. Open symbols are for an isotropic emission in the co-moving frame, filled symbols 
are for an anisotropic case: radiation emitted within 47r/17 sr of the radial flow direction. Note that the 
second peak (T p ~ 3 s) shows a stronger increase in spectral hardness in the latter case and that in both 
cases the maximum spectral hardness occurs ~ 0.5 — 1 s before the intensity peak, (c) The first peak as seen 
in each BATSE channel, (d) The second peak in the same bands. Fluxes in (c) and (d) are normalized to 
the peak value in that channel. Legends in (c) and (d) give the peak time, the duration (the sum of the rise 
and fall times, as defined in text), and the time-asymmetry ratio of each peak, in this order, in each BATSE 
channel. 
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Fig. 6. — BATSE window time histories of a burst with F = 400 and different degrees of anisotropy of 
radiation emission in the co-moving frame. Solid curves: isotropic emission (two cones, each of Att/2 sr 
solid angle, corresponding to a cap of angular opening T _1 in the laboratory frame). Dotted and dashed 
curves: anisotropic emission (47r/5 sr and 47r/17 sr co-moving frame solid angles, corresponding to 0.5 F -1 
and 0.25 r _1 angular opening caps, respectively). Note that, as the co-moving frame anisotropy increases, 
the light-curves becomes more symmetric. Legends indicate the rise and fall times and the asymmetry ratio, 
as defined in text, in this order. Upper graph: photon fluxes, lower graph: energetic fluxes. Light-curves 
have been normalized to their peak fluxes, binned in 256 ms, and aligned at their peaks. 
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Fig. 7.— Tq = 400 multi-peak burst (in BATSE window) generated by a variable magnetic field (inset of 
upper graph). Upper graph: isotropic emission. Middle and lower graph: anisotropic emission. The solid 
angle around the direction of motion in which the emission is confined is indicated in each graph. Note that 
as the degree of anisotropy increases, the peaks become more distinct. Photon and energy fluxes have been 
normalized at the maximum value reached during the burst, aligned at their peak times, and binned in 64 
ms. Numbers indicate the expected beginning of each pulse. 



